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Abstract. The measurement and characterization of the lensing of the cosmic microwave 
background (CMB) is a key goal of the current and next generation of CMB experiments. 
We perform a case study of a three-channel balloon-borne CMB experiment observing the sky 
at (Z,6)=(250°,— 38°) and attaining a sensitivity of 5.25 ^K— arcmin with 8' angular resolution 
at 150 GHz, in order to assess whether the effect of polarized Galactic dust is expected to be a 
significant contaminant to the lensing signal reconstructed using the EB quadratic estimator. 
We find that for our assumed dust model, polarization fractions as low as a few percent may 
lead to a significant dust bias to the lensing convergence power spectrum. We investigate a 
parametric component separation method, proposed by Stompor et al. (2009), for mitigating 
the effect of this dust bias, based on a 'profile likelihood' technique for estimating the dust 
spectral index. We find a dust contrast regime in which the accuracy of the profile likelihood 
spectral index estimate breaks down, and in which external information on the dust frequency 
scaling is needed. We propose a criterion for putting a requirement on the accuracy with 
which the dust spectral index must be estimated or constrained, and demonstrate that if this 
requirement is met, then the dust bias can be removed. 
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1 Introduction 

The measurement and characterizaton of the weak lensing of the cosmic microwave back- 
ground (CMB) by the large-scale structure distribution is a promising and active field of 
research in observational cosmology (for a review of the physics of CMB weak lensing, see 
[1]). Measurements of this signal can break fundamental degeneracies that afflict the cos- 
mological interpretation of measurements of the CMB power spectrum [2] as well as help to 
improve the constraints on the cosmological parameters [3, 4]. As a result, a number of on- 
going and planned experiments are targeting the weak lensing signal as one of their primary 
science goals. 

The CMB weak lensing signal was first detected by [5] who cross correlated data from 
the Wilkinson Microwave Anisotropy Probe (WMAP) satellite with a tracer of large-scale 
structure in the form of the NRAO VLA sky survey. Since then observational progress has 
been rapid: the Atacama Cosmology Telescope (ACT) collaboration made the first detec- 
tion of weak lensing signal using CMB data alone [6] and the South Pole Telescope (SPT) 
collaboration have followed with a detection at higher significance [7] as well as detecting 
the correlation of the weak lensing 'convergence' and large-scale structure tracers from the 
Wide- field Infrared Survey Explorer and Spitzer/IRAC [8]. First applications of the weak 
lensing signal measurements have been to provide corroborating evidence for the cosmological 
constant from CMB data alone [9] and to improve constraints on the dark energy equation 
of state [7]. In the future, improved cosmological constraints are expected from the full SPT 
and ACT surveys, and especially from Planck which is poised to significantly advance lensing 
studies [10, 11]. 

The most widely adopted technique for extracting the lensing signal, also implemented 
in this investigation, is the 'quadratic estimator' suggested by [12] which is a near optimal 
framework for reconstructing the lensing field using a quadratic combination of an appropri- 
ately filtered CMB temperature and/or polarization map and its gradient. The subsequent 
estimation of the power spectrum of the lensing field is met by a variety of real world com- 
plications. These include intrinsic biases of the quadratic estimator which must be modeled 
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and subtracted [11, 13], and any effect that breaks the statistical isotropy of data including 
anisotropic noise [14] and the effect of instrumental systematics [15, 16]. 

Astrophysical foregrounds represent another source of contamination to the lensing sig- 
nal for which there have been several simulation studies assessing their possible impact and 
suggesting mitigation strategies. [17] performed simulations of the kinetic Sunyaev-Zeldovich 
(SZ) effect and pointed out that this non-Gaussian effect correlated with the lensing signal, 
with a spectrum indistinguishable from CMB anisotropies, could potentially bias the lensing 
reconstruction. They proposed masking CMB maps around the location of clusters detected 
via the thermal SZ effect as a potential strategy for mitigating kSZ biases, a procedure subse- 
quently adopted by the SPT team [7] . [18] made calculation of the lensing biases induced by 
polarized extragalactic radio point sources, concluding that the effect of these source could 
be mitigated via source detection/masking, and that these sources ought not pose a signif- 
icant challenge to a future CMB polarization satellite. [10] studied the impact of diffuse 
and extragalactic foregrounds on a Planck-like simulations and proposed a multi-frequency 
masking and 'in-painting' technique for measuring the lensing signal. 

To date there has been no specific study on the possible impact of Galactic polar- 
ized dust emission on the detection of the lensing signal. We believe that it is important to 
make such an assessment given the fact that several ongoing ground-based and balloon-borne 
CMB polarization experiments including ACTPol [19], SPTPol [20], EBEX [21] and POLAR- 
BEAR [22] are targeting the lensing signal in the near future, using measurements of CMB 
polarization in the frequency range 90-410 GHz. In particular we will focus on performing 
a case study of an EBEX-like experimental configuration which, owing to observing con- 
straints on balloon-borne telescopes flying from Antarctica, will most likely survey a region 
of sky where Galactic foregrounds in the form of polarized dust emission may be significant 
compared to the lensing signal. We aim to develop methods for numerically investigating 
the extent to which this experiment can expect to be challenged by diffuse polarized dust 
emission, and to investigate possible dust mitigation strategies. 

This study is set out as follows: Section 2 describes the data analysis techniques that 
we have implemented; Section 3 describes the CMB and dust polarization simulations that 
we have performed; Section 4 describes our findings and results, and Section 5 contains our 
conclusions. 

2 Data analysis techniques 

In this section we review the formalism of the two analysis techniques that we have imple- 
mented and investigated in our lensing estimation pipeline. These are the quadratic estima- 
tor [12] for estimating the CMB lensing signal, and the multi-frequency analysis technique 
of [23] for estimating the foreground spectral index and performing separation of the CMB 
component. 

2.1 CMB lensing and the quadratic estimator 

CMB photons from the last scattering surface are deflected according to 





(Q + iU)(n) = (Q±iU) n + V^(n) 
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where 4>(n) is the projected potential, and n is the direction on the sky [24]. We will adopt 
the notation convention that lensed quantities are denoted with a tilde and unlensed ones 
without. In the Born approximation, in which the deflection angles are assumed to be small 
enough to carry out the projection kernel integration along the line of sight, eft is related to 
the three-dimensional gravitational potential, ip, by 

#(n) = -2 j dD </, (On.O), 

(2.2) 

where D and D* are the comoving angular diameter distances to the lens and the CMB last 
scattering surface respectively; <p is related to the 'convergence' k as V 2 (p = [24]. 

To linear order in eft, the changes in the Fourier moments of the temperature and polar- 
ization fields due to lensing are [25] 

Sf(l) = J -0^T(l')W(l',L), (2.3) 

5E(l) = J ^2^(l / )cos2 Wl - J B(l , )sin2^i]w(l , ,L), 

8B(l) = J -^-[B(l')cos2 Wl + E(l') S m2<p vl ]w(l',L), 

where (fy\ = ipy — <pi, L = 1 — 1', and 

^(1,L) = -[1-I#(L). (2.4) 

From these equations we note that the effect of lensing is to couple the gradient of the 
primordial CMB 1' modes to that of the observed 1 modes. Furthermore, starting from zero 
primordial S-modes, B(Y) = 0, lensing generates I?-mode anisotropies in the observed map. 

Quadratic estimators for the convergence: We will assume a CMB map with homo- 
geneous white noise and a Gaussian beam smoothing. The power spectrum of the detector 
noise is [26] 

Cf = a 2 pix O pix , (2.5) 

where o" p i x is the RMS noise per pixel and O p i x is the solid angle subtended by each pixel. The 
observed CMB temperature and polarization fields, X G [T,E,B], and their power spectra, 
, are 

X? hs = li e- 1 ^^ + Nf, (2.6) 
C?>° hs = Cf e~ l2 °l + Cf^, 

where is the Fourier mode of the detector noise, and relates to the full-width half- 
maximum (FWHM) of the telescope beam, ^fwhm, via #fwhm = o"6 V 7 8 In 2. 

To extract the lensing information in the observed CMB map we will use the quadratic 
estimator formalism [12, 27, 28] in the context of the convergence estimators [29, 30]. These 
estimators are uniquely determined by the requirement that each estimator be unbiased 
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(k XY (n)) = «(n) over an ensemble average of the CMB temperature and polarization fields 
X and Y, and the variance of the estimator be minimal, 

(k X Y k* v XY ) = (2vr) 2 5 D {1 - l')(Cf + JVp xy ). (2.7) 

In real space the convergence estimators are given by [30] 

G XY (n) - ilXy ^pbi| e 2mj e 2 6 . ( 2 - 8 ) 



Wv(n) = /y^o 4fe^ \e-* 1 (2.9) 



d 2 l Yj obs f e 2 ^ 1 

where is the azimuthal angle of the wavevector 1; the two phase factors in braces are 
applied when Y = E, B respectively, and is unity when Y = T. Also C XY = C XE for 
Y = B. The construction of these fields incorporate the deconvolution of the beam from the 

1/2 2 

map, hence the beam factors e~2 °6 appearing on both fields. 

Given the two filtered fields in Eq. (2.8) and Eq. (2.9), the convergence estimators are 
then given by 



aXY n 

kf Y = l — il ■ / d 2 n Re [G X Y(n)Wy(n)] e~ ll il . (2.10) 

The normalization coefficients, A XY , are related to the noise power spectrum, N?' X , of the 
estimators k XY (h) by Nf' XY = l 2 A XY /A, and are calculated as 



1 f d\ (l-h)C XY f XY 



A XY I 2 J (2vr) 2 c X ' obs C Y ' obs 
\ sin2Av9 J 

with 1 = li + 1 2 , A</j = ip h - tpi 2 , and pf h Y l2 ) = f x £ <fa , where [28] 



(2.11) 



/il,i 2 = (1-11)^ + (1-12)^, (2.12) 

= (l'h) Cgc^Atp + Q-h) C%, 

/iIJ 2 = (l-li)Q?sm2A^, 

fhf, = [(Mi)Qf + (l-l 2 )Qf]cos2A^, 

/ 1 f,f 2 = (l-li)Qfsin2A^. 

Our code for estimating the convergence using the quadratic estimator formalism is a direct 
implementation of the above equations, Eq. (2.5)~(2.12). 

2.2 Foreground spectral index estimation and cleaning 

To perform foreground spectral index estimation and cleaning, we use an implementation of 
the parametric component separation algorithm proposed by [23] and tested in the context of 
forecasts for inflationary i?-mode detection in [31, 32]. In this framework, the multi-frequency 
data vector d p is modeled at each pixel p of the map as 

dp — ApSp -\- Tip, (2.13) 
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where A p = A p ((3) is an A^f rC q x A^ com p 'mixing' matrix with. A^ S p CC free parameters, to be 
estimated, s p is a vector of N comp component amplitudes also to be estimated, and n p is the 
noise. A likelihood for the data is given by 

- 21n£(s,/3) = CONST + (d — As) t N~ 1 (d - As), (2.14) 

where iV _1 is the noise covariance of the map. A key result of [23] is that under the as- 
sumption that the spectral index is the same for each pixel, the maximum-likelihood values 
of the components s can be found by first locating the values of the spectral parameters (3 
that minimize the value of the 'profile likelihood' 

- 21n£ spec (/3) = CONST - (A t N~ 1 d) t (A 1 N^ 1 A)" 1 (A^^d) , (2.15) 

which is an expression independent of s. Once the maximum likelihood spectral parameters, 
(3, have been determined, their values are substituted into the generalized least squares 
solution of Eq. (2.13), given by 

s = (A t N~ 1 A)~ 1 A*JV _1 d, (2.16) 

N s = (A'AT _1 A) _1 , (2.17) 

to obtain the estimated component amplitudes, s, their noise covariance Ng, pixel by pixel. 

Insight into the component separation error can be gained using the ^-matrix formalism 
derived in [31] where, 

Z0) = (A\P)N- l A0)) ~* A\p)N~ l A (fa) . (2.18) 

In our case, this will be a 2 x 2 matrix corresponding to the two components, CMB and dust. 
In the limit of /? = /3q then Z is the identity matrix. For the case where the dust spectral 
index is mis-estimated, then the off-diagonal terms of Z quantify the fraction of the original 
component that remains unsubtracted, since 

f-comp 

§i= Zii0)*j(Po), ( 2 -19) 

3 

where s are the estimated components. 

3 CMB and polarized dust simulations 

We simulate CMB polarization and diffuse polarized dust emission on a 13° x 13° patch of 
sky located at (RA, Dec) = (75°, -44.5°), corresponding to (/,6)=(250°,-38°) in Galactic 
coordinates. This area of sky is accessible from observing sites both in Antarctica and 
Chile, has been surveyed by several past ground-based and balloon-borne CMB experiments 
including Boomerang [33], QUAD [34], ACBAR [35], QUIET [36] and SPT [37], as shown in 
Figure 1. It is close to the area that will be observed by the balloon-borne CMB experiment 
EBEX, owing to the fact that is aligned with the anti-sun direction in mid December, soon 
after the start of the Antarctic long duration balloon flight launch window. Given that the 
anti-sun direction moves by 1 degree in RA per day in the direction of the Galactic plane, 
then the sky patch we assume can be considered representative of an early launch scenario 
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Figure 1. Dust intensity at 150 GHz (thresholded to 25 /^Krj). Black outlines show the borders of 
previous CMB intensity and polarization surveys, while the white outline shows the 13° x 13° patch 
that we consider in this study. 



when the impact of Galactic foregrounds will be at their lowest. Following [32], we assume a 
three band experimental configuration with channels at 150, 250 and 410 GHz observing to 
depths of 5.25, 14.0 and 140 /xKcmb - arcmin respectively, each with an angular resolution 
of 8'. An advantage of this selection of bands is to minimize the possible impact polarized 
synchrotron, assumed to be negligible in this study, which may affect ground-based CMB 
polarization experiments observing at 90 GHz [31]. 

Our Galactic polarized dust model is the same as the one first described in [31]: dust 
intensity is given by the model of [38] extrapolated to 410 GHz. Then, to simulate po- 
larized emission, polarization angles are set on large angular scales using the WMAP dust 
template [39], while on smaller scales, extra Gaussian power is added using the prescription 
of [40]. The polarization fraction, p, is assumed to be spatially constant, and we investigate 
three cases of 3.6, 5, and 10%, intended to bracket the average high Galactic latitude dust 
polarization detected in the WMAP W band [41], and possible higher dust polarization frac- 
tions observed by ARCHEOPS at 353 GHz [42]. The dust is scaled from the 410 GHz band 
to the lower frequency bands assuming a greybody frequency scaling 

Aiust oc j- -, (3.1) 

exp$-l 

with T = 18K and /3 = 1.65, with the dust temperature and spectral index both assumed to 
be uniform across the patch. The resulting dust polarization simulation at 150 GHz is shown 
in Figure 2. 

For our CMB simulations, we produced two sets of 100 realizations-lensed and unlensed- 
with 0.76' pixel size, assuming the WMAP 7- year best-fit cosmological parameter values [43], 
and with our fiducial CMB polarization power spectra calculated using CAMB [44]. 




1.57 4.79 8.00 11.22 14.43 17.65 -0.42 -0.25 -0.08 0.08 0.25 0.42 -0.57 -0.34 -0.11 0.11 0.34 0.57 



Figure 2. Dust simulation (Stokes /, Q, and U) at 150 GHz, with a polarization fraction p = 0.036, 
on a 13° x 13° patch centred on (Z,6)=(250° -38°). 




-9.53 -5.72 -1.91 1.91 5.72 9.53 -7.24 -4.34 -1.45 1.45 4.34 7.24 -8.66 -5.19 -1.73 1.73 5.19 8.66 



Figure 3. From left to right, we show our CMB+dust+noise simulations at 150, 250 and 410 GHz 
for a dust polarization fraction p = 0.1. The 150 GHz channel is CMB dominated, while the 410 GHz 
channel is dust dominated. The 250 GHz channel provides information about the dust spectral index. 

The first set of maps-the lensed CMB realizations-were obtained starting from the 
unlensed fiducial power spectrum C^ Y , from which Gaussian realizations of the CMB po- 
larization were generated, which were then lensed by remapping the pixels by the deflection 
field. The deflection field is in turn is derived from a Gaussian realization of the projected 
potential power spectrum Cf^; we neglect the effect of the integrated Sachs- Wolfe effect in- 
duced correlation C^. We have chosen our pixel size to be small compared to the RMS 
of the deflection angles (~ 2') so that errors due to interpolation back onto the regular grid 
after remapping are small. We have checked that the E and .B-mode power spectra of these 
simulated lensed maps reproduces the lensed power spectra obtained from CAMB to within 
a few percent accuracy for the .E-mode spectrum and to within five percent accuracy for 
the .B-mode spectrum. While this is less accurate than the all-sky lensing simulations now 
performed by several groups using various interpolation schemes [45] [46] [47] [4], we believe 
that our flat-sky simulations are sufficiently accurate for our dust foreground study. 

The second set of maps-the unlensed CMB realizations-were obtained from Gaussian re- 
alizations of the lensed fiducial power spectrum, Cf Y . These maps have the same power spec- 
trum as the lensed CMB realizations, but have none of the lensing-induced non-Gaussianity. 
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Figure 4. Validation of the EB quadratic estimator. The solid line is the power spectrum of the input 
convergence map C^ K , and the error bars show the estimated convergence power spectra averaged 
over 100 simulations. The dash curves are the lensing noise spectrum, N?' , for the two sets of 
survey parameters assumed. 



Since the power spectrum of the convergence reconstructed on unlensed CMB maps is same 
as the lensing noise power spectrum predicted analytically from Eq. (2.11), these maps have 
been used for checking the accuracy and implementation details of the convergence and power 
spectrum estimators, as well as the testing effect of mask apodization. 

Finally the CMB maps are scaled to antenna temperature units in the three bands at 
150, 250 and 410 GHz, smoothed with an 8' beam, and uncorrelated Gaussian white noise is 
added to each pixel. The Stokes Q parameter of an example simulation is shown in Figure 3. 

4 Results 

This section describes our results in which we calculate the level of lensing bias that is 
expected from our dust polarization model. Finding a dust bias to be present, we apply 
the parametric multi- frequency foreground cleaning technique [23] described in Section 2.2, 
and establish a requirement for the accuracy with which the dust spectral index must be 
constrained in order to guarantee dust cleaning. 

4.1 Foreground-free case 

Before reporting the effect of the dust, we first demonstrate our reconstruction of the con- 
vergence from the EB quadratic estimator, using modes in the range 100 < £ < 3000, under 
the most idealized foreground-free case. Figure 4 shows the convergence power spectrum for 
our simulated 150 GHz channel with 5.25 fiK— arcmin sensitivity and 8' angular resolution, 
as well as for a survey with 4^tK— arcmin sensitivity and 1.4' angular resolution similar to the 
planned 'ACTPol Deep' survey of [19]. The power spectrum estimates shown are the average 
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Figure 5. The lensed CMB .E-mode (left panel) and -B-mode (right panel) power spectrum is 
compared to the power spectrum of the dust (jp = 0.1) and instrumental noise at 150 GHz. For 
comparison the fiducial E and .B-mode power spectra is also shown (cyan). 



over 100 simulations of C^ K — N^ XY , while the corresponding error bar is given by 

ACr = l , =A , (4.1) 

where we have applied a binning scheme Al = 98 (thirty bins between I = 40 and I = 3000). 
The power spectrum estimates shown in Figure 4 represent an end to end validation of our 
CMB simulations, EB quadratic estimator, and power spectrum estimation pipeline. 

4.2 Dust polarization bias at 150 GHz 

To first assess the size of the dust contamination on the patch we are considering, we estimated 
the power spectrum of the simulated dust at 150 GHz and compared it to E and .B-mode 
signal and noise power spectra, as shown in Figure 5 for the example of p = 0.1. For our dust 
model and choice of patch, the E and .B-mode power spectra of the dust approximately follow 
a powerlaw given by Cf nst = (A x p) 2 i" , where p is the polarization fraction, A ~ 120/^K 
and /? ~ —3.5. Our previous study [32] has shown that polarized dust at this level of power 
must be modeled and subtracted in order to derive unbiased estimated of the inflationary B- 
mode spectrum, a cosmological signal which is accessible in the £ < 200 range of the .B-mode 
power spectrum. The main question we seek to address in this study is whether this level of 
anisotropy power of foreground contamination is large enough to also bias the estimates of 
the lensing signal. 

We have calculated the power spectrum of the convergence field reconstructed with the 
EB quadratic estimator using the dust contaminated 150 GHz channel, and show our results 
in Figure 6 for three different polarization fractions of p = [0.01,0.036,0.1]. We find that 
if the dust contamination is ignored during the lensing estimation, then a 'dust noise bias' 
dominates over the lensing power spectrum estimates for both the p = 0.036 and p = 0.1 
case, but can be ignored for the p = 0.01 case. 
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Figure 6. Effect of polarized dust at 150 GHz for surveys with sensitivity A p = 5.25/xi^— arcmin. 
The upper, middle and lower panels show cases with dust polarization fractions of p = 0.001, 0.036 
and 0.1 respectively. We find that for our dust model and choice of patch, diffuse polarized dust is 
expected to be a significant source of bias to if the dust polarization fraction is greater than 
around 1 percent. 



Although our demonstration of the dust bias will be dependent on the choice of patch 
we have assumed, and on the details of our polarized dust model and its power spectrum, 
we nonetheless conclude that diffuse polarized dust may in principle be a source of bias for 
future sub-orbital CMB surveys aiming at lensing estimates using the EB estimator, and 
that methods for foreground debiasing must therefore be developed. 

For our dust model and choice of patch we can now roughly estimate that the require- 
ment on the 'dust suppression factor' (the factor by which the polarized dust must is reduced 
at 150 GHz) is approximately (pcrit/p), where we have estimated that the 'critical polariza- 
tion fraction' p cr it) below which the effect of dust contamination becomes small compared to 
the noise level in the convergence power spectrum estimates, is approximately p cr a = 0.01. 
This requirement on the dust suppression factor, (p C rit/p)> wm m turn set the requirement 
with which the dust spectral index must be estimated, which we will calculate in the next 
section. 

4.3 Debiasing the effect of dust via foreground cleaning 

Having established that our polarized dust model is leading to a bias in the lensing power 
spectrum, our main aims are 1) To test the multi-frequency foreground cleaning technique 
described in Section 2.2, in which the dust spectral index is first estimated from the data 
using the profile likelihood Eq. (2.15), before applying a linear least-squares component sep- 
aration of the CMB and dust using Eq. (2.16) and propagation of the noise covariance using 
Eq. (2.17), and 2) To put requirements on the accuracy with which the dust spectral index 
must be estimated in order to guarantee dust cleaning. 

For the homogeneous noise case that we have simulated, the final least squares combina- 
tion of the data is a linear combination of the channel maps. However, the profile likelihood is 
a non-linear method involving quadratic combinations of the data, and so we have therefore 
numerically investigated its accuracy in determining the dust spectral index by carrying out 
100 spectral index estimation and component separation simulations, each with a different 
CMB and noise realization. 

Table 1 summarizes our results for our dust spectral index estimation simulations, for 
three cases p = [0.036,0.05,0.1]. First, we find that the greater the dust signal, the more 
accurately the dust spectral index can be estimated from the data. This is qualitatively 
consistent with the analysis of the profile likelihood of [48] in which increased dust contrast 
is shown to improve the accuracy of the dust spectral index estimation. We have also found 
that the profile likelihood returns an unbiased estimate for the dust spectral (3 when consid- 
ering the average value obtained over the ensemble of CMB+noise realizations that we have 
analysed, but we also find that the profile likelihood width for any single noise underestimates 
the dust spectral index uncertainty. This breakdown of the profile likelihood appears to be 
related to the relatively low contrast regime of the polarization data we have simulated: If 
the dust intensity data is assumed to have the same frequency scaling as the polarization 
data, then the dust spectral index is first estimated with much higher accuracy, and the 
profile likelihood width and Monte Carlo average of the spectral index estimates agree well. 
This implies that any application of the profile likelihood to data should be accompanied 
by Monte Carlo simulations, given a Galactic foreground model, in order to more accurately 
quantify the spectral index uncertainty. 

For each of the 100 spectral index estimates, we calculate the Z-matrix defined in 
Eq. (2.18); We will focus on the ^cMB,Dust matrix element, since this will quantify the frac- 
tion of dust that remains unsubtracted from the CMB component after the least squares 
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(Anput = 1-65) 



p 




Ap 


0.036 


1.69 


0.43 (0.04) 


0.05 


1.67 


0.22 (0.03) 


0.1 


1.66 


0.06 (0.02) 



Table 1. Mean and RMS of the dust spectral index estimate /3, estimated from 100 component sep- 
aration simulations varying CMB and noise realizations, for the three dust polarization fraction cases 
assumed; The values in parentheses are the mean value of the profile likelihood width. We find that 
spectral index uncertainty for any single noise realisation (shown in parentheses) is underestimated 
by the profile likelihood in the dust contrast regime of our polarization simulations. 



p = 0.036 



Input: 


CMB 


Dust 


Output: 






CMB 


1 


-0.06 ±0.79 


Dust 





1.04 ±0.41 


p = 0.05 


Input: 


CMB 


Dust 


Output: 






CMB 


1 


-0.01 ±0.41 


Dust 





1.01 ±0.21 


p = 0.1 


Input: 


CMB 


Dust 


Output: 






CMB 


1 


0.00 ±0.10 


Dust 





1.00 ±0.05 



Table 2. Z-matrices, Eq. (2.18), averaged over 100 component separation simulations varying CMB 
and noise realizations, for three different dust polarization cases. The -ZcMB.Dust matrix element 
quantifies the fraction of dust at 150 GHz that is mixed into the CMB after component separation, 
and so can be thought of as a 'dust suppression factor'. We find that in the p = 0.1 case, the dust 
spectral index is estimated with enough accuracy to guarantee dust cleaning (| (ZcMB,Dust) | <C 1), 
while for the p = 0.36 case dust cleaning is not guaranteed. 

component separation. Table 2 gives the mean and RMS of Z-matrices for the three polar- 
ization fraction cases p = [0.036,0.05,0.1]. The basic conclusion is that while in the case 
of p = 0.1 the dust spectral index is always estimated accurately enough to guarantee dust 
cleaning ((ZcMB,Dust) = 0.00 ± 0.10), in the case of p = 0.036 the dust is no longer bright 
enough to allow sufficiently accurate spectral index estimation, and for many simulations, 
there is an amplification of the dust contamination (|ZcMB,Dust| > !)■ 
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Figure 7. Convergence power spectrum of the estimated CMB component. The upper and 
middle panels shows the average convergence power spectrum for realisations with spectral index 
/3 > [1.81,1.76,1.71] and /3 < [1.49,1.54,1.59], respectively, for the three dust polarization fraction 
cases of p = [0.036,0.05,0.1]. The lower panel shows the average convergence power spectrum for 
those CMB components whose spectral index meets the accuracy requirement given in Eq. (4.2). 
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It therefore becomes important to be able to put a requirement on the accuracy, A/3, 
with which the dust spectral index must be estimated in order to guarantee foreground 
cleaning. This requirement can then either guide the application of our parametric com- 
ponent separation technique, or inform the usage of possible external spectral index prior 
information. For our case study we suggest the criterion 



P 



(4-2) 

0o 



where we conservatively assume p CT n = 0.01; More generally, the factor (pcrit/p) can be 
thought of as the desired dust suppression factor for a given patch of sky. Qualitatively, 
there is therefore a requirement for greater spectral index precision for larger values of p. 



Quantitatively, we have numerically calculated - 



= 0.56 from the function Z(B), and so 

00 



dz 

for the three cases p = [0.036,0.05,0.1] we require A/3 < [0.16,0.11,0.06] respectively. We 
have checked that enforcement of this spectral index accuracy requirement leads to satisfac- 
tory dust cleaning. Figure 7 shows the average convergence power spectrum of the estimated 
CMB component for those realizations whose spectral index estimate meet this accuracy 
requirement (for instance those with /3 = 1.65 ± 0.06 for the p = 0.1 case), for which we find 
satisfactorily debiased convergence power spectrum estimates after component separation. 
Conversely, the realizations with estimated dust spectral indices falling outside the required 
range show a biased convergence power spectrum after component separation. 



5 Conclusions 



Several ongoing and planned CMB polarization experiments are aiming to measure and 
characterise the lensing of the cosmic microwave background, in order to improve constraints 
on the parameters of the cosmological model. Within this context we have made the first 
specific study of the possible effect of diffuse polarized dust emission on the accuracy of the 
reconstruction of the lensing convergence signal. Our particular focus has been on performing 
a case study of a three channel balloon-borne CMB experiment covering the frequency range 
150-410 GHz. Our numerical investigation is based on a dust polarization simulation and 
a flat-sky implementation of the Hu and Okamoto quadratic estimator. We found that for 
the sky patch under consideration, which is near to the region of sky that will be targeted 
by the EBEX experiment, and for plausible dust polarization fractions in the range 3.6- 
10%, the anisotropy of the diffuse dust polarization will be large enough at 150 GHz to 
bias the reconstruction of the convergence. Thus a multi-frequency experimental approach 
is imperative, and appropriate analysis methods must be developed for debiasing the effect 
of polarized dust. 

In order to mitigate the effect of the dust and to debias the convergence power spec- 
trum, we demonstrated a multi-frequency component separation technique in which the dust 
spectral index is first estimated from the data using a profile likelihood technique, before ap- 
plying a least squares component separation. We found evidence for a dust contrast regime 
in which the accuracy of profile likelihood breaks down, both underestimating the spectral 
index uncertainty as well as providing spectral index estimates that are insufficiently accurate 
to guarantee dust cleaning. This highlights the possible need for external constraints on the 
frequency scaling of the polarized dust which can then be used as a prior, to stabilize the dust 
spectral index estimates to the accuracy required for sufficient dust cleaning. We proposed a 
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criterion, Eq. (4.2), which sets a requirement for the accuracy with which the spectral index 
of the foregrounds must be estimated in order to guarantee a given dust suppression factor. 
We then demonstrated that satisfactory dust cleaning was achieved for the cases in which 
the estimated spectral index met this requirement. 
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